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Abstract 

In this work we explore the effects of human mobihty on the dispersion of a vector borne disease. 
We combine an already presented stochastic model for dengue with a simple representation of the 
daily motion of humans on a schematic city of 20x20 blocks with 100 inhabitants in each block. 
The pattern of motion of the individuals is described in terms of complex networks in which links 
connect different blocks and the link length distribution is in accordance with recent findings on 
human mobility. It is shown that human mobility can turn out to be the main driving force of the 
disease dispersal. 

PACS numbers: 87.10.Mn, 87.23. Ge, 05.10.Gg 
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I. INTRODUCTION 



Dengue fever is a vector borne disease produced by a flavivirus of the family flaviviridae 
The main vectors of dengue are Aedes aegypti and Aedes albopictus. 

Dengue fever epidemics have been known for more then two centuries as benign fever 
epidemics appearing at intervals of one to several decades. Recently the number of countries 
affected by this epidemic has increased and severe forms have become more frequent. 

The research aimed at producing dengue models for public policy use began with Newton 
and Reiter [2^ who introduced a minimal model for dengue in the form of a set of Ordinary 
Differential Equations (ODE) for the human population disaggregated in Susceptible, Ex- 
posed, Infected and Recovered compartments. The mosquitoes population was not modeled 
in this early work. A different starting point was taken by Focks et al. js, 4] that began by 
describing mosquitoes populations in a computer framework named Dynamic Table Model 
where later the human population (as well as the disease) was introduced j^. 

Newton and Reiter's model (NR) favors economy of resources and mathematical acces- 
sibility, in contrast. Pock's model emphasize realism, these models represent in Dengue two 
contrasting compromises in the standard trade off in modeling. A third starting point has 
been recently added. Otero and Solari (OS) developed a dengue model [6|] which includes 
the evolution of the mosquitoes population [T], |8| and is spatially explicit. This last model is 
somewhat in between Pock's and NR as it is formulated as a state-dependent Poisson model 
with exponentially distributed times. 

ODE models have received most of the attention. Some of the works explore: variability of 
vector population 9j, human population 10|, the effects of hyp othetical vertical transmission 
of Dengue in vectors Q, seasonality [J age structure R as well as incomplete gamma 
distributions for the incubation and infectious times P]. Comparison with real epidemics 
has shown that there is a need to consider the spatial heterogeneity as well [ISj . 

In a previous work we have developed a dengue model which includes the evolution 
of the mosquitoes population and is spatially explicit. In that work the spatial spread of the 
infection was driven by the flights of the mosquitoes that gave rise to a diffusion process. In it 
we analyzed the evolution of dengue infection in a city of 20X20 blocks with 100 individuals 
in each one, this population was fixed throughout the calculation and no mobility of the 
individuals was allowed. 
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As such, in that model, the spatial evolution of the dengue infection was only driven by 
the flight of mosquitoes as the mobility of humans was not included. It is usually recognized 
that human mobility is not only necessary to be included in human infection spread models, 
but that it might be the main source of the dynamics behind spatiotemporal phenomena 
on geographic scales (i.e the spread of infection from city to city due to people flying long 
distances by plane). It is thus very important to address the problem of the mobility of 
humans and incorporate it into the models to be able to make more reliable predictions, 
and then, to be able to propose effective public policies against the dispersal of a known or 
emerging disease. 



Including the mobility of the human population in a model is not an easy task given the 
complexity of human behavior. The flrst problem to address is the technical and ethical 
difficulties that arise when trying to get information about the mobility of humans. There 
are many databases from which this data could be inferred, such as the ones associated 
to cellular phone networks, credit cards, hotel reservations, flight reservation databases, 
etc. But as almost all of them are private, most researchers do not have access to them. 
Moreover, even if we did have them, mixing this diversifled information together to get a 
human mobility model is a hard task by itself. Aside from this particular difficulties, there 
is an intrinsic bias on the databases if we are going to use them for diseases spread, because 



its reasonable to think t 
awareness of a disease 17 



lat human behavior will change, or adapt in presence of social 



26| . and the inferences made on this databases can not take that 



into account. Moreover there can be a social bias because not everyone may use credit cards, 
go to hotels, etc. 

Whether it is necessary to have a detailed information on the movements of each indi- 
vidual to build up a model, or if it is only needed a coarse grain statistics of the mobility as 
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a whole is an open question wich still has to be answered 

Several works tackle the issue of the correct description of the human beings mobility, 
relaying on different methods and databases 281-135 1. 



As most works focusing on t 



human transmitted diseases 



lis topic analyze the effect of the human mobility in 
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numan- 
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36l438l| , but not in vector borne ones 27 
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in the present work we show a variation of our previous dengue model which not only includes 
the flight of the mosquitoes but also the mobility of the humans beings. We then show the 
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key differences in the results between both our models and conclude on the actual impact 
of the human behavior on dengue. 

In section II we describe the characteristics of the epidemiological model we use in this 
work encompassing the dynamics of the virus for humans, mosquitoes and the dispersal 
dynamics for each. In section III we present the results of our numerical investigations 
which include the analysis of the size and time evolution of the epidemics, the morphological 
properties of the patterns of spatial distributions of the infections for different mobility 
patterns and for different densities of mosquitoes. Finally conclusions are drawn in section 
IV 



II. THE EPIDEMIOLOGICAL MODEL 

There are four ingredients in this model, the epidemiological dynamics of the infected 
mosquitoes, the epidemiological dynamics of the infected humans and the mobility pattern 
of the individuals and mosquitoes. Each of this elements will be discussed in what follows 



A. Mosquitoes 

The dengue virus does not make any effect to the vector, as such, Aedes aegypti pop- 
ulations are independent of the presence of the virus. In the present model mosquitoes 
populations are produced by the Aedes aegypti model Q| with spatial resolution of one 
block using climatic data tuned to Buenos Aires, a temperate city where dengue circulated 
in the summer season 2008-2009 4l|]. The urbanistic unit of the city is the block, approx- 
imately a square of (100m x 100m). Because of the temperate climate the houses are not 
open as it is often the case in tropical areas. Mosquitoes develop in the center of the block 
which often presents vegetation and communicates the houses of the block. The model then 
assumes that mosquitoes belong to the block and not to the houses and they blood-feed 
with equal probability in any human resident in the block. Aedes aegypti is assumed to 
disperse seeking for places to lay eggs. The mosquitoes population, number of bites per day, 
dispersal flights and adult mortality information per block is obtained from the mosquitoes 
model 
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The time step of the model has been fixed at one day. 

The virus enters the mosquito when it bites a viremic human with a probabihty Phm.{j) 
depending of the day j in the infectious cycle of the human bitten. The cycle continues with 
the reproduction of the virus within the mosquito (extrinsic period) that lasts days. After 
this reproduction period the mosquito becomes infectious and transmits the virus when it 
bites with a probability pmh- The mosquito follows a cycle Susceptible, Exposed, Infected 
(SEI) and does not recover. Eventually mosquitoes die with a daily mortality of 0.09 3]- 
The adult female mosquitoes population as produced by the Aedes aegypti simulation is 
then split into susceptible, stages of exposed and one infective compartment according to 
their interaction with the viremic human population and the number of days elapsed since 
acquiring the virus. 

The mosquito population of each block is not fixed, but instead mosquitoes move around 
in terms of a simple diffusion process. 



B. Humans 

The evolution of the disease in one individual human, /i, evolves as follows: 
Day 

The virus is transmitted by the bite 
d = do — )■ of an infected mosquito. 

The human enters the exposed stage 

d = do + TE^h) — )► The human h enters the infective stage 

The human is infective and 
d = do + TE^h) + j — > transmits the virus to a biting mosquito 

with probability Phm{j) 



The human enters the recovered stage. 

d > do + teQi) + Tj -> 

No longer transmits dengue 
In the table above TE{h) stands for the intrinsic incubation time, and Xj is the viremic time 
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of each individual. The cycle in the human being is then of the form Susceptible, Exposed, 
Infected, Recovered (SEIR). Each human has it's own value of te which is assigned according 



to the Nishiura's experimental distribution 42 1. 



As mentioned above, our analysis of the time evolution of dengue fever is performed on 
a schematic city in which the basic unit is the block and in each block a human population 
of 100 individuals is placed. An exposed human(index case) is introduced near the center 
of the grid, on January 1**. 

The human population of each block is not fixed in the present work. In order to describe 
the patterns of mobility of the humans we have adopted the following schematic model. 50% 
of the population of each block is randomly selected to be mobile, while the other 50% is 
considered to remain in its original block during the whole analysis. Each mobile individual is 
assumed to stay 2/3 of the day in its original block, while the other 1/3 of the day she/he will 
stay in a randomly assigned block according to the corresponding distribution. Each of the 
mobile individuals is assigned in each case a fixed destination to which it returns everyday. 
At the end of the day individuals return to their original block. This random assignment is 
performed according to certain rules that will characterize the mobility pattern. Following 
recent works on human mobility, referred in the introduction, we require that the movement 



le lengths 
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of each individual a) should be highly predictable [43j] and b) the distribution of t 
of the displacements of the human should follow a truncated Levy distribution 
reads: 

P(r) oc (r + ro)~'^exp(— r/zt). (1) 

Being P(r) the probability of a human traveling a distance, where tq, /3 and k are parameters 
that characterize the distribution. In this work we have used the parameters described in 
Table E 

TABLE I: Levy-Flight distribution parameters. 

ro(m) P K (m) 

1 100 1.65 1500 

2 100 2. 1500 

3 100 3. 1500 

4 100 4. 1500 
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Such a pattern of mobility of the humans is accomphshed by building a network with 50 
links starting in each block. The length of the link is distributed according to the proposed 
length distribution and the final block is chosen at random from those which can be reached 
by the link. Each link is assigned to a mobile human at the start of the simulation. The 
distribution of jumps sizes for each type of underlying network is shown on Fig[T] 

In order to have reference mobility networks we have also analyzed the case in which 
the endpoints of the links are completely random. We have also investigated the case in 
which only one individual per block performs a random jump while the rest of the mobile 
individuals in the block visit only their neighboring blocks {MoveAOO). 
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FIG. 1: Distribution of jump length for different mobility patterns. Random(Full circles), 
Move400(Empty circles), Levy-Flights with (3 = 1.65(Full squares), f3 = 2(Full triangles), 
(3 = 3(Crosses), f3 = 4(Full diamonds). As comparison, with asterisk is shown the Levy-Flight 



with parameters from the Barabasi distribution 



Once the set of parameters is fixed, different underlying networks are generated and a set 
of evolutions (typically a couple of thousand events) is performed for such arrangements. 



The networks 



The networks built according to the above mentioned prescription can be analyzed in 
order to unveil their "small world properties". We have found it interesting to study the 
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geodesic path i.e. the average minimum path between all the cells. 



^E^^.- (2) 



n{n 

with dij the minimum path between cells i and j, and n the total number of nodes. The 
minimum path is defined as the minimum number of links that are to be traversed in order 
to travel from the original block to the destiny block. Therefore, we see that it is a simple 
average over all the possible pairs of blocks in the system of the minimum distance between 
each pair. 

Another interesting magnitude to explore the characteristics of a network is the so called 
Clusterization. One of the usual definitions of this magnitude is: 
Given a node i with ki nearest neighbors we define q as 

2 



(number of links between nearest neighbors) (3) 



^ = ^E-^ (4) 

i 

In Table [III we show the results of such a calculation. 

TABLE II: Clusterization and Mean shortest path for the human mobility networks. 

Network I C 

Random 1.779 0.220 

Levy(/3 = 1.65) 2.006 0.286 

Levy(/3 = 2) 2.112 0.303 

Levy(/3 = 3) 2.532 0.349 

Levy(/3 = 4) 3.092 0.380 

Only 400 move 3.884 0.014 



We can see from table I that the mean shortest path attains a minimum for the completely 
random pattern of links and grows as this pattern is replaced by the ones generated by the 
levy fiights. We see that the broader the Levy-Flight is the larger is the average minimum 
path, as expected. 
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III. NUMERICAL CALCULATIONS 



We have implemented the above described method (The algorithm without human mobil- 
ity has been fully described in ^6\) and have performed extensive calculations for different 
initial conditions. 

The conditions are : mosquitoes breeding sites density, different realizations of the un- 
derlying mobility networks and different seasonal conditions. 

The number of breeding sites per block explored in this calculations are 50, 100,200,300 
and 400. Larger number of breeding sites are considered to be too unrealistic for the system 
we have in mind i.e. the city of Buenos Aires. Moreover larger number of breeding sites do 
not add new information to our calculations. 

Another relevant condition that we have explored is the underlying mobility network. 
For each value of the breeding sites density, evolutions with different underlying networks 
were performed. 

Finally we have considered two different seasonal conditions, on Fig. [2] we show the 
population(top) and temperature(bottom) profiles. In the season-less situation, the temper- 
ature remains fixed all along the evolution at 23°C. In this case the mosquito population 
remains basically constant all along the evolution. In such a case the size of the epidemics 
is determined by the dynamics of the infection itself subject to the above mentioned bound- 
ary conditions. If we adopt the average temperature time-distribution of BA (see figure 
for details) the population of mosquitoes is a strongly time dependent one, the size of the 
epidemics might then be severely constrained. 



In what follows we will focus on certain properties of the epidemic system that are relevant 
for the understanding of the characteristics of the time evolution. In first place we study the 
morphology of the evolving spatial structure of the epidemics. Then we study the final size 
and time span of the epidemics. Then we include the results of the analysis of a system in 
which the temperature is kept fixed at 23 degrees Celsius. Finally we study a new magnitude 
that we name the power of the epidemics. 
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FIG. 2: Top: Normalized mosquitoes population for constant temperature(Dashed) and seasonal 
variation(full). Bottom: Temperature profile for constant and seasonal variation. Day is set at 
1^* of January. 

A. Morphology of the spatial structure of the epidemics 

It is expected that human mobility increases the size and speed of the epidemics. 

This happens because each jump (shortcut) when executed by a virus carrying individual 
may induce the contagion of mosquitoes at the destination block and then generates a new 
dispersal center for the illness. 

In Fig [3] we show the density of recovered individuals at three relevant times for the case 
in which the dispersal of dengue is driven by the diffusion of mosquitoes only. It is seen that 
the population of recovered individuals displays a symmetrical pattern as expected from a 
simple diffusion process. 

On the contrary as seen in Fig H] the pattern for the case in which the human jumps 
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FIG. 3: Spatial distribution of the number of recovered individuals for three times namely 25 
days, 81 days and 173 days for the case of a dengue epidemic driven only by the dispersal of 
mosquitoes. It can be seen that the pattern corresponds to a diffusive process and the evolution is 
highly symmetric. 

follow a Levy-Flight distribution is quite heterogeneous and it can be clearly seen that there 
is more than one dispersal center. 

This observations can be made more quantitative if we calculate the radial correlation 
function defined as the probability of finding at least one infected (recovered) individual 
(calculated at the time at which all individuals have returned home) at a block such that it 
can be reached by a jump of length r from the place at which the initial infected individual 
one was located (which in this case is a block close to the center of the city). 

In what follows we show the result of calculating the radial correlation function for three 
typical cases namely for the mosquito driven evolution, for the case Levy-Flight with /3 = 3 
and the case in which only one of the mobile individuals performs random jumps. 

From Figs. I3|^|71 is clearly seen that in the presence of human mobility nearly all of the 
city can be reached by the epidemic in short times. On the one hand the case of mosquitoes 
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FIG. 4: For the same times as in Fig. [31 the spatial distribution of recovered individuals when 
their mobility is modeled by a Truncated Levy-Flight with (/3) = 3. It can be seen that, at variance 
with Fig. [3] the pattern is asymmetrical, moreover the evolution is faster and it can be seen that 
at t = 81 days a second focus is present. 

only as driving force the correlation function displays patterns expected for a traveling 
wavefront. In the other cases the wavefront breaks early in the evolution and the correlation 
function is different from zero almost everywhere after a few days. 

As we have seen in Fig. 3) and 4) the structure of the spatial density of, say, humans in 
state R is highly symmetric and compact for the case without human mobility. As human 
mobility (of the kind considered in this work) is incorporated both the symmetry and the 
compactness are lost. In order to explore this behavior in a more quantitative way we define 
cluster of recovered individuals in the following way. Given a block i we will call it an 
occupied block if at least one member of its original population is in the recovered state. A 
cluster (of size larger that one) is a set of occupied blocks in which all constituents have at 
least a nearest or second nearest neighbor which belongs to the cluster. Then the block i 
will belong to the cluster if the following relation is satisfied: 
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Distance from index case 

FIG. 5: Radial correlation function as defined in the text for the case in which the spatial evolution 
of the epidemics is driven by the dispersal of mosquitoes only. Full squares (red)t = 25days, 
full triangles (green)t = 57days, crosses (blue)t = 89days, diamond (violet)t = 121days open 
circles (black) t = 153days, open squares (gray) t = 189days. The pattern of variation is quite 
regular as expected for a diffusive case. 

i E C ■<=^ 3 j 6 C / i is neighbor of j (5) 
We define the mass of a cluster as the number of recovered individuals in the cluster. 

The results are displayed in the Fig [HI 

It can be seen that for the case of simple diffusion all of the mass is concentrated in the 
biggest cluster . On the other hand for the random case there is a time (around 50 days into 
the epidemic) at which only about 65% of the mass is in the biggest cluster . This is due to 
the emergence of secondary foci generated by infective humans who perform long jumps. 

B. Sizes and time span of the epidemics 

One of the main observables in this kind of problems is the final size of the epidemics. In 
what follows we show Figure a comparison of the final size of the epidemics in terms of the 
boxplots corresponding to different values of the breeding sites density for different patterns 
of human mobility. As described above we have two limiting situations a) the case in which 
the moving humans perform jumps with completely random destinations and /) the case in 
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FIG. 6: Radial correlation function as defined in the text for the case in which the humans move 
with the Levy-Flight (/? = 3) distribution. Full squares (red)t = 25days, full triangles (green)t = 
bldays, crosses (blue)t = 89days, diamond (violet)t = 121days open circles (black) t = 153days, 
open squares (gray) t = 189days. 

which the dispersal of the epidemics is only due to the diffusion of mosquitoes. In between 
we have the patterns related to the length of the jump given by a truncated Levy-Flight 
characterized by the different set of parameters shown in table [T] i.e. /3 = 1.65,2,3 and 4 
(keep in mind the smaller /3 the closer to the random case). The corresponding results are 
displayed in panels b — e. 

Finally in Fig. [TD]we show (left panel) the case in which only one of the mobile humans in 
each block performs a jump to a random destination while the others move to nearest neigh- 
bors. For the sake of completeness we show on the right panel the boxplot corresponding to 
the mosquito only driven evolution. 

It is immediate to see that the effect of human mobility for all the cases is to increase 
the final size of the epidemics with respect to the case in which the mosquitoes diffusion 
is the only driving force. Moreover in the case of completely random mobility and for the 
Levy-Flight with /3 = 1.65 we get that for the highest BS density proposed in this work, the 
epidemics spreads over the whole population. 

Figure [H] shows the duration distribution of the epidemics as a function of the different 
patterns of human mobility and for two constant breeding site densities of 200 BS/ha (Top) 
and 400 BS/ha (Bottom) respectively. For 200 BS (Top) all boxplots present a similar 
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FIG. 7: Radial correlation function as defined in the text for the case in which only one human per 
block moves at random, and the rest move to neighboring blocks. Full squares (red)t = 25days, 
full triangles (green)t = 57days, crosses (blue)t = 89days, diamond (violet)t = 121days open 
circles (black) t = IbSdays, open squares (gray) t = 189days. 

spread of data but a slight tendency of increase of the median from left to right (from the 
case without mobility to the complete random pattern). Instead, for 400 BS the tendency 
is to decrease from left to right. (This behavior will be properly discussed in terms of the 
Power of the epidemics, see below) 

C. Behavior of the model at constant temperature 

Figure [12] shows the duration distribution of the epidemics as a function of the different 
patterns of human mobility for a constant breeding site density of 400 BS/ha and a constant 
temperature of 23 degrees Celsius. We see that the maximum duration of the epidemic takes 
place for the case of no human mobility and it shortens as the patterns of mobility of the 
individuals tend to the completely random one. It is also interesting to note that in this 
case the epidemics involves the whole population. 
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FIG. 8: Ratio between the size of the biggest cluster (in terms of the number of recovered humans) 
and the total number of recovered. Diffusive (Red full squares), random (Green full triangles), 
Levy-Flight(/3 = 1.65) (blue crosses), Levy-Flight (/3 = 4) (violet full diamonds), Move400 (Black 
asterisks) . 

D. Power of the epidemic 

We define the mean power of the epidemic as the ratio between the median of the final 
size of the epidemic and the median of the duration of the epidemics. 

Pm = - (6) 

T 

Figure [I3] shows the mean power for three conditions: (a) 400 BS and constant temper- 
ature of 23 degrees, (6) 400 BS and seasonal variation of temperature and (c) 200 BS and 
seasonal variation of temperature. The Pm grows with broader jump length distributions 
and with higher BS densities. For the case of 400 BS is higher for constant temperature 
than for seasonal variation of temperature. If we compare 400 BS and 200 BS (for seasonal 
variation of temperature) we see the same increasing tendency of Pm with human mobility 
but this value is higher for 400 BS than for 200 BS. (For the case of 400 BS and constant T 
the increase of the power of the epidemic is a consequence of the reduction of its duration 
as the pattern of human mobility approaches the fully random case. At constant T all of 
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FIG. 9: Boxplots graphs of the epidemic size for different patterns of human mobihty: (a) random, 
(b) levy(/3=1.65), (c) \evy{P=2), (d) levy(/9=3), (e) levy(/3=4) and (f) move400. 

the population gets infected. On the other hand, when the Temperature is not fixed it 
severely constrains the mosquitoes population and then the increase in is mainly due to 
the increase of the size of the infected population as the time span of the epidemics is only 
mildly dependent on the driving force of the dispersal. 



IV. CONCLUSIONS 

In this work we have explored the effect of human mobility on the dynamics of a vector 
borne infection. Wc have added this characteristic of human behavior on an already tested 
model of dengue dispersal when the dynamics is driven by mosquitoes alone. We have 
analyzed the case of a schematic city of 20x20 blocks with 100 individuals per block. 

We have considered two temperature profiles, on the one hand a simple constant tem- 
perature one and a realistic time distribution corresponding to the city of Buenos Aires, 
Argentina. 
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FIG. 10: Boxplots graphs of the epidemic size for different patterns of human mobihty: (Left) 
move400 and (Right) no human mobihty. 



Another variable in our analysis has been the number of breeding cites in the city, we 
have considered 50,100, 200, 300 and 400 breeding cites per block. 

Human mobility has been described by superimposing diverse kinds of networks in which 
links represent the daily movement of humans. The distribution of lengths of this links 
are derived from recent studies on human motion and in particular, taking into account 
the finding that human behavior is highly predictable. We have also considered reference 
patterns i.e purely random motion and random motion of a single human per block. 

We have explored different observables like size and duration of the outbreaks, and com- 
plementary the morphological characteristics of the pattern of recovered individuals. 

We have found that human mobility strongly enhaces the infection dispersal. Even for 
the case in which just one individual per block can perform a long jump. This effect can be 
traced to the fact that when the disease dispersal is driven by mosquitoes alone we have a 
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FIG. 12: Comparison of the duration of epidemics for the different patterns of human mobility at 
a constant temperature of 23 degrees Celsius (400 breeding sites / ha) . 

single focus that expands due to diffusive kind of diynamics. When human mobility is taken 
into account, multiple loci appear as the time evolution is followed 

Human mobility increases the size and the speed of propagation of the outbreaks. This 
feature can be captured by the magnitude " Power of the epidemics" defined as the quotient of 
the size of the epidemics divided by the its time span. This magnitude displays a monotonous 
increase as the Mean Length Path of the network describing the daily human mobility pattern 
decreases. 

This findings indicate that human mobility might turn out to be the main driving force 
in the epidemics dynamics. 

Both in the case of fixed temperature and seasonal variational one, human motion gives 
rise to faster and more widespread epidemics. 

Finally this findings indicate that, when considering meassures to fight epidemics dispersal 



20 



300 




Average Minimum Path (1) 



FIG. 13: Comparison of the mean power of the epidemics for the different patterns of human 
mobihty (from left to right: Random, Levyl.65, Levy2, LevyS, Levy4, Move400), characterized 
by the average minimum path for three conditions: (a) 400 BS and constant temperature of 23 
degrees and (6) 400 BS and (c) 200 BS for seasonal variation of temperature. 

human motion should be one of the top concerns. We are presently exploring this issue. 
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